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^ ■ Recent measurements of J /if) production in e+e~ colliders pose a chal- 

\ lenge to the NRQCD factorization theorem for quarkonium production, 

ff^ i Discrepancies between leading order calculations of color-octet contribu- 

' tions and the momentum distribution of J/ -0 observed by Belle and BaBar 

are resolved by resumming large perturbative and nonperturbative cor- 
rections that are enhanced near the kinematic endpoint. The large cross 
sections for J/ij^cc and double quarkonium production remain poorly un- 
derstood. Nonperturbative effects in fixed-target hadroproduction of open 
' charm are also discussed. Large asymmetries in the production of charm 

mesons and baryons probe nonperturbative corrections to the QCD factor- 
^ ^ ization theorem. A power correction called heavy-quark recombination can 

5^ ' economically explain these asymmetries with a few universal parameters. 
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1. Introduction 



In the last couple of years there have been a number of interesting exper- 
imental results in the production of heavy particles, including measurements 
of J/ijj production in e+e" collisions at the T(4S') resonance ^ HI 01 and 
charm meson production at the Tevatron This talk focuses on how these 
results impact theoretical understanding of heavy particle production. The 
spectrum of J/ip observed in e^e~ colliders disagrees with leading order 
calculations based on Non-Relativistic QCD (NRQCD) factorization the- 
orems 1^1 . Better agreement is obtained in calculations which resum the 
large nonperturbative and perturbative corrections that arise near the kine- 
matic endpoint [6.. However, large cross sections for J j^cc and exclusive 
double quarkonium production remain poorly understood. I also discuss 
heavy quark production. Calculations of charm and bottom production 
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which resum large logarithms of p± / mq provide a consistent description of 
the production of these particles at the Tevatron 7, 8 , 9 . Finally, I explain 
how the large asymmetries observed in fixed-target hadroproduction exper- 
iments can be explained by a power correction to the QCD factorization 
theorem called heavy-quark recombination lOi illj . 



The current theoretical framework for understanding the production of 
heavy quarkonia is NRQCD 0. NRQCD solves important theoretical and 
phenomenological problems in quarkonium theory. Color-singlet model cal- 
culations of Xc decay suffer from infrared divergences Jj^. NRQCD provides 
a generalized factorization theorem that includes nonperturbative correc- 
tions to the color-singlet model, including color-octet decay and production 
mechanisms. Infrared divergences are factored into nonperturbative ma- 
trix elements, so calculations of inclusive production and decay rates are 
infrared safe JH]- Color-octet production mechanisms are necessary for un- 
derstanding the production of J/tp at large transverse momentum, p±, at 
the Tevatron Jl]. Convincing evidence for color-octet mechanisms has re- 
cently been seen in 77 collisions, where color-singlet mechanisms underesti- 
mate the cross section by an order of magnitude while calculations including 
color-octet mechanisms describe the data well jl5j . 

However, there are many unsolved problems in quarkonium physics |16j . 
Perhaps the most puzzling is the polarization of J/V' and ip' , which is pre- 
dicted to be transverse at very large p± in hadron colliders jiTj. The theo- 
retical prediction is consistent with the data at intermediate p± but at the 
largest p± measured the J/tp and ip' are observed to be slightly longitudi- 
nally polarized. 

This talk focuses on new puzzles arising from recent measurements of 
J/ip production in e"'"e~ colliders ^ |H The leading color-octet con- 
tribution to this process was expected to dramatically enhance the cross 
section for maximally energetic J/t/j as well as modify their angular distri- 
bution 18 . If cos 9 is the angle of the J/ip relative to the axis defined by 
the e~^e~ beams and z = Ej/^/Ey,^, then the differential cross section is 



The function A{z) tends to -1 as z ^ 1 for color-singlet production. The 
leading color-octet diagram contributes only at z = 1 and gives ^(1) ~ 1. 
The total color-singlet cross section is predicted to be 0.4-0.9 pb ^^1; while 
the total cross section from the leading color-octet mechanism is expected 
to be ~ 1 pb. A substantial rise in the cross section near the kinematic 
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endpoint accompanied by a change in angular distribution was predicted to 
be a robust signal of the color-octet mechanism ^18^. 

Experimental data does not agree with these expectations One 
problem is that a sharp rise in the cross section near the kinematic endpoint 
is not observed. On the other hand, the cross section is larger than predicted 
by the color-singlet model and A{z) ~ 1 for 0.7 < z < 1. Another puzzle is 
the production of open charm with J/ip. Belle observes that a large fraction 
of J/ tp are produced with an additional cc P : 



a{e+e -^J/jjcc) _ +0.15 



0.59;[]:}^ ±0.12. (2) 



(Preliminary results described by P. Pakhlov in this conference suggest that 
the ratio is even larger [20] •) Leading order color-singlet model calculations 
predict the ratio to be « 0.2 jT^] and a large color-octet contribution will 
make the ratio even smaller. 

The resolution of the first problem lies in a careful analysis of the per- 
turbative and nonperturbative corrections that appear near the kinematic 
endpoint of quarkonia production [HJI^. The J/^p production cross section 
in NRQCD is 

^(e+e- ^J/^P + X)=Y: ^(e+e- - cc[n] + , (3) 

where {On'^) are NRQCD matrix elements and da{e^e~ cc[n] + X) are 
perturbatively calculable short-distance cross sections. The label n denotes 
the angular momentum and color quantum numbers of the cc. The NRQCD 

scaling rules show that {On^) scales as some power of v, where v is the 
typical velocity of the heavy quarks inside the bound state. 

The leading Feynman diagrams which contribute to color-octet produc- 
tion give da{e'^e~ —>■ cc[n] + X)/dz proportional to 6{1 — z). This is the first 
in an infinite series of terms that are singular at z = 1. There are higher 
order nonperturbative corrections that scale as — z)"" j21j as well as 

perturbative corrections that go like a"ln™'(l — z)/{l — z), m < 2n — 1. For 
charmonium, ~ ~ 0.3 so for z > 0.7 perturbation theory and the 
NRQCD V expansion both break down. 

Near the kinematic endpoint, the final state consists of two kinds of 
quanta: energetic collinear particles with light-like momenta in the jet 
against which the J/'0 is recoiling and particles which are soft as viewed 
from the rest frame of the J/ip- NRQCD breaks down because the theory 
does not explicitly include these degrees of freedom. The problem can be 
fixed by using the Soft-Collinear Effective Theory (SCET) |^ which ex- 
plicitly includes both collinear and soft degrees of freedom. By combining 
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Fig. 1. The sum of the resummed color- octet and leading order color-singlet con- 
tributions are plotted as the upper line. The lower line is the leading order color- 
singlet contribution only, and the data are from the BaBar collaboration (left) 
and Belle collaboration 'Ji] (right). 



SCET and NRQCD one finds that in the endpoint region Eq. Q is replaced 
with the following factorization theorem: |2j 

'^ = at^ [diS^-\i)J{i-z). (4) 

Here cjg is a short distance cross section which is perturbatively calculable. 
The shape function S'""! (^) is a universal nonperturbative distribution that 
resums the large nonperturbative corrections. — z) is a calculable func- 
tion that describes the propagation of the collinear particles in the energetic 
gluon jet. Large perturbative corrections can be resummed by solving the 
SCET renormalization group equations for Cq"^, 5["](^) and — z). 

Comparison of this calculation with data from BaBar ^ and Belle jSj is 
shown in Fig. ^ ^ . The number of events is plotted as a function of the J/ if^ 
momentum, p^. The lower line in both plots is the leading color-singlet con- 
tribution which falls well below data. The upper line includes the resummed 
color-octet cross section. The leading order calculation of the color-octet 
contribution (not shown) gives an integrated cross section comparable to 
that of the resummed calculation shown here, but the entire cross section 
is located in the very last bin of . This also does not agree with data. 
Fig. n demonstrates that when large perturbative and nonperturbative cor- 
rections are included the momentum spectrum of the J/ij: produced by the 
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color-octet mechanism is significantly broadened. The agreement with data 
is good in the endpoint region where the calculation is most reliable. The 
calculation is not predictive because the shape function is fitted to available 
data. However, the moments of the shape function can be estimated using 
the NRQCD scaling rules and the shape function used in Ref. ^ satisfies 
these constraints. The universality of the shape function can be tested by 
applying the methods of Ref. [S] to other J/il) production processes. 

While resolving the discrepancy between leading order color-octet cal- 
culations and the observed distribution, the calculation does not help 
explain the large cross section for J /tpcc observed by the Belle collabora- 
tion. The Belle collaboration also observes a large cross section for exclusive 
J/ip + r]c and J/ip + Xc [H]- An NRQCD analysis of exclusive cross sections 
appears in Ref. [22]. Because the helicity conservation rules for exclusive 
QCD processes suppress the leading QCD contribution, the purely QED 
contributions are surprisingly large (~ 20%). The leading relativistic cor- 
rections give large corrections which unfortunately are difficult to estimate 
reliably. For instance, Ref. |131 quotes (T{J/'ip + r]c) = 5.5'tl%^ fb, with 
the uncertainty dominated by relativistic corrections. For comparison Belle 
measures a{J/ip + ric) xBilrjc ^ 4 charged particles] = 33lg±9 fb. Clearly a 
mechanism for enhancing the J /ipcc and double charmonium cross sections 
is needed. Proposals for nonperturbative mechanisms responsible for this 
enhancement appear in Ref. j24j . 

3. Open Charm and Bottom Production 

The QCD factorization theorem states that the production cross section 
for a heavy particle H containing a heavy quark Q is j^Hl 

da[AB ^ = E f^/A ® fj/B ® da[ij ^ QQX] ® Dq^h + - • (5) 

Here /j/^ is a parton distribution function, Dq-,h is a fragmentation func- 
tion, da[ij — > QQX] is a short-distance cross section and the ellipsis repre- 
sents corrections which are suppressed by Aqcd/w-q or Aqcd/p±- 

Perturbative aspects of Eq. © are tested by measurements of charm 
and bottom production at collider experiments. Experimental reviews of 
heavy quark production at LEP, HERA and the Tevatron are described in 
the talks of A. Sciaba [^Hl) B. Olivier |27, and K. Rinnert respectively. 
At the Tevatron, discrepancies between NLO calculations of bottom produc- 
tion and experimental cross sections have been known for a long time j29j . 
Recently CDF has extended measurements to include charm as well as bot- 
tom |1] . Resummation of large logarithms of p±_ / mq is needed to obtain 
better agreement with both charm and bottom cross sections [TjlHlini. It 
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is also important to treat fragmentation functions carefully as the total 
cross section is sensitive to the fragmentation function used [HI- Though 
the existing calculations differ in how finite mass corrections are handled 
they are consistent numerically and agree with data within the theoretical 
uncertainties estimated by varying factorization and renormalization scales. 

Nonperturbative power corrections to Eq. @ are probed in lower en- 
ergy fixed-target experiments. Production asymmetries are an incisive test 
of these corrections. At leading order in perturbation theory, charm par- 
ticles and antiparticles are produced symmetrically because the partonic 
processes gg cc and qq — > cc produce charm and anticharm symmetri- 
cally and Dc^H = ^c-^Tl ^° charge conjugation invar iance. At next- 
to-leading order, the asymmetry, a[H] = {(t[H] — a[H\) / {a[H] + a[H\), is 
only a few percent [30] . Fixed-target hadroproduction OH ESI IS^ and 
photoproduction |351 136j experiments observe much larger asymmetries. In 
hadroproduction the asymmetries are known as the "Leading Particle Ef- 
fect". Cross sections for charm particles sharing a valence quark with the 
beam hadron are enhanced in the forward direction of the beam. Hadropro- 
duction asymmetries can be quite large. For example, in the most forward 
region measured in 7r~A^ collisions, the ratio of D~ to is ~ 6. 

Charm asymmetries are conventionally explained by nonperturbative 
models of hadronization |37j . These models suffer from a lack of predictive 
power, since they depend on a number of nonperturbative functions, such 
as the distribution of spectator quarks in hadron remnants. The most com- 
monly used model is the Lund string fragmentation model j39j which can be 
implemented using PYTHIA [10]. The PYTHIA Monte Carlo with default 
parameters rarely predicts the asymmetries correctly |31| and in the case of 
Ac asymmetries in -kN collisions [SH and collisions |36j gets the sign of 
the asymmetry wrong. 

A novel mechanism for generating charm hadron asymmetries called 
heavy-quark recombination has recently been introduced in Ref. [10|. Simi- 
lar mechanisms for production of light hadrons were considered in Ref. [U . 
An important difference between the production of heavy hadrons and light 
hadrons is that in the former case heavy quark symmetry |42jj can be used 
to simplify the structure of nonperturbative factors appearing in the cal- 
culation. Heavy-quark recombination is an 0(Aqcd/"Ic) suppressed power 
correction to the factorization theorem of Eq. Q. In this process, a light 
anti-quark, g, from the incident hadron participates in a hard-scattering 
process which produces a c and c quark. Following the hard scattering the 
q and the c recombine to form a D meson. A similar mechanism in which 
a light quark recombines with the c quark is the dominant recombination 
contribution to charm baryon production Heavy-quark recombination 
differs from previous nonperturbative models in that the asymmetry is gen- 
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Fig. 2. Asymmetry as a function of xp for mesons (top left) and D*^ mesons 
(top right) produced in 500 GeV nN collisions '33! and mesons produced in 
600 GeV Yi~ N collisions 'JUl (bottom). Theory curves are explained in text. 



erated in the short-distance process so cross sections are calculable up to an 
overall normalization that is set by a few universal parameters. The short- 
distance cross section is strongly peaked in the forward direction of the light 
quark or antiquark, naturally leading to an asymmetric cross section. 

Heavy-quark recombination accounts for the D meson asymmetries ob- 
served in photoproduction and hadroproduction experiments. Fig. |21 shows 
asymmetries for and D*~^ mesons produced in 500 GeV vrA^ collisions 
\6'6\ and asymmetries for Ds mesons produced in 600 GeV Ti~N collisions 
j34j . There are four universal parameters in the theory. The theory curves 
in Fig. [21 correspond to fits with one, two and three of these parameters. (A 
four parameter fit yields identical results as the three parameter fit.) The 
heavy-quark recombination mechanism correctly describes asymmetries for 
different types of D mesons in experiments with different beams with a 
minimal set of universal parameters. The heavy-quark recombination also 
correctly describes Ac asymmetries in both vrA^ and pN collisions jllj . 



8 



REFERENCES 



[I] K. Abe et al. [Belle Collaboration], Phys. Rev. Lett. 89, 142001 (2002) 

[2] B. Aubert et al. [BABAR Collaboration], Phys. Rev. Lett. 87, 162002 (2001). 
[3] K. Abe et al. [BELLE Collaboration], Phys. Rev. Lett. 88, 052001 (2002). 
[4] D. Acosta et al. [CDF Collaboration], hep-ex/0307080, 

[5] W. E. Caswell and G. P. Lepage, Phys. Lett. B 167, 437 (1986); G. T. Bodwin, 
E. Braaten and G. P. Lepage, Phys. Rev. D 51, 1125 (1995); ibid. D 55, 5853 
(1997); M. E. Luke, A. V. Manohar and L Z. Rothstein, Phys. Rev. D 61, 
074025 (2000). 

[6] S. Fleming, A. K. Leibovich and T. Mehen, Phys. Rev. D 68, 094011 (2003). 
[7] J. Binnewies, B. A. Kniehl and G. Kramer, Phys. Rev. D 58, 034016 (1998). 
[8] M. Cacciari and P. Nason, Phys. Rev. Lett. 89, 122003 (2002). 
[9] M. Cacciari and P. Nason, JHEP 0309, 006 (2003). 

[10] E. Braaten, Y. Jia and T. Mehen, Phys. Rev. Lett. 89, 122002 (2002); Phys. 
Rev. D 66, 034003 (2002). ihid., 014003 (2002). 

[II] E. Braaten, M. Kusunoki, Y. Jia and T. Mehen, "hep-ph/0304280l 

[12] R. Barbieri, R. Gatto and E. Remiddi, Phys. Lett. B 61, 465 (1976); Phys. 
Lett. B 106, 497 (1981); R. Barbieri, M. Caffo, R. Gatto and E. Remiddi, Nucl. 
Phys. B 192, 61 (1981). 

[13] G. T. Bodwin, E. Braaten and G. P. Lepage, Phys. Rev. D 46, 1914 (1992). 

[14] E. Braaten and S. Fleming, Phys. Rev. Lett. 74, 3327 (1995); P. L. Cho and 
A. K. Leibovich, Phys. Rev. D 53, 150 (1996); Phys. Rev. D 53, 6203 (1996). 

[15] M. Klasen, B. A. Kniehl, L. N. Mihaila and M. Steinhauser, Phys. Rev. Lett. 
89, 032001 (2002). 

[16] G. T. Bodwin, |hep-ph/0212203| 

[17] P. L. Cho and M. B. Wise, Phys. Lett. B 346, 129 (1995); A. K. Leibovich, 
Phys. Rev. D 56, 4412 (1997); M. Beneke and M. Kramer, Phys. Rev. D 55, 
5269 (1997); E. Braaten, B. A. Kniehl and J. Lee, Phys. Rev. D 62, 094005 
(2000). 

[18] E. Braaten and Y. Q. Chen, Phys. Rev. Lett. 76, 730 (1996). 

[19] P. L. Cho and A. K. Leibovich, Phys. Rev. D 54, 6690 (1996); F. Yuan, 
C. F. Qiao and K. T. Chao, Phys. Rev. D 56, 321 (1997). S. Back, P. Ko, 
J. Lee and H. S. Song, J. Korean Phys. Soc. 33, 97 (1998). G. A. Schuler, Eur. 
Phys. J. C 8, 273 (1999). J. H. Kuhn and H. Schneider, Phys. Rev. D 24, 2996 
(1981); Z. Phys. C 11, 263 (1981); V. M. Driesen, J. H. Kuhn and E. Mirkes, 
Phys. Rev. D 49, 3197 (1994). 

[20] P. Pakhlov, proceedings of this conference. 

[21] M. Beneke, I. Z. Rothstein and M. B. Wise, Phys. Lett. B 408, 373 (1997). 

[22] C. W. Bauer, S. Fleming and M. E. Luke, Phys. Rev. D 63, 014006 (2001). 
C. W. Bauer, S. Fleming, D. Pirjol and L W. Stewart, Phys. Rev. D 63, 
114020 (2001). C. W. Bauer and L W. Stewart, Phys. Lett. B 516, 134 (2001). 



9 



C. W. Bauer, D. Pirjol and I. W. Stewart, Phys. Rev. Lett. 87, 201806 (2001); 

Phys. Rev. D 65, 054022 (2002); Phys. Rev. D 66, 054005 (2002). 
[23] E. Braaten and J. Lee, Phys. Rev. D 67, 054007 (2003). 
[24] B. L. loffe and D. E. Kharzeev, hep-ph/0306062, A. B. Kaidalov, JETP Lett. 

77, 349 (2003) [Pisma Zh. Eksp. Tcor. Fiz. 77, 417 (2003)]. 
[25] J. C. Colhns, D. E. Soper and G. Sterman, Nucl. Phys. B 263, 37 (1986). 
[26] A. Sciaba, proceedings of this conference. 
[27] B. Ohvier, proceedings of this conference. 
[28] K. Rinnert, proceedings of this conference. 

[29] D. Acosta et al. [CDF CoUaboration] , Phys. Rev. D 65, 052005 (2002). 

[30] P. Nason, S. Dawson and R. K. Elhs, Nucl. Phys. B 327, 49 (1989) ; 
W. Beenakker, H. Kuijf, W. L. van Neerven and J. Smith, Phys. Rev. D 40, 
54 (1989); W. Beenakker, W. L. van Neerven, R. Meng, G. A. Schuler and 
J. Smith, Nucl. Phys. B 351, 507 (1991). 

[31] M. Adamovich et al. [WA82 CoUaboration], Phys. Lett. B 305, 402 (1993); 
G. A. Alves et al. [E769 Collaboration], Phys. Rev. Lett. 72, 812 (1994); ibid. 
77, 2392 (1996); E. M. Aitala et al. [E791 Collaboration], Phys. Lett. B 411, 
230 (1997); M. Adamovich et al. [BEATRICE Collaboration], Nucl. Phys. B 
495, 3 (1997); M. I. Adamovich et al. [WA89 Collaboration], Eur. Phys. J. C 
8, 593 (1999); F. G. Garcia et al. [SELEX Collaboration], Phys. Lett. B 528, 
49 (2002). 

[32] E. M. Aitala et al. [E791 Collaboration], Phys. Lett. B 495, 42 (2000). 
[33] E. M. Aitala et al. [E791 Collaboration], Phys. Lett. B 539, 218 (2002); ibid. 
371, 157 (1996); 

[34] M. Kaya et al. [SELEX CoUaboration], Phys. Lett. B 558, 34 (2003), 
[35] J. .C. Anjos et. al. [Tagged Photon Spectrometer Collaboration], Phys. Rev. 
Lett. 62, 513 (1989); M.P. Alvarez et. al. [NA14/2 Collaboration], Z. Phys. C 
60, 53 (1993); P. L. Frabetti et. al. [E687 CoUaboration], Phys. Lett. B 370, 
222 (1996); J. C. Anjos and E. Cuautle [FOCUS Collaboration], AIP Conf. 
Proc. 531, 172 (2000). 
[36] J. M. Link et. al. [FOCUS Collaboration], hep-ex/0311022, 

[37] R. Vogt and S. J. Brodsky, Nucl. Phys. B 478, 311 (1996); E. Norrbin and 
T. Sjostrand, Phys. Lett. B 442, 407 (1998); A. K. Likhoded and S. R. Sla- 
bospitsky, Phys. Atom. Nucl. 62, 693 (1999) [Yad. Fiz. 62, 742 (1999)]. 
O. I. Piskounova, hep-ph/0202005, R. C. Hwa, Phys. Rev. D 51, 85 (1995). 
R. Rapp and E. V. Shuryak, Phys. Rev. D 67, 074036 (2003). T. Tashiro, 
S. Nakariki, H. Noda and K. Kinoshita, Eur. Phys. J. C 24, 573 (2002). 
E. Cuautle, G. Herrera and J. Magnin, Eur. Phys. J. C 2, 473 (1998). 

[38] G. Herrera and J. Magnin, Eur. Phys. J. C 2, 477 (1998) 

[39] H.-U. Bengtsson and T. Sjostrand, Comput. Phys. Commun. 46, 43 (1987). 

[40] T. Sjostrand, L. Lonnblad and S. Mrenna, hep-ph/0108264, 

[41] K. P. Das and R. C. Hwa, Phys. Lett. B 68, 459 (1977) [Erratum-ibid. 73B, 
504 (1978)]; E. L. Berger, T. Gottschalk and D. W. Sivers, Phys. Rev. D 23, 
99 (1981). 



10 



[42] N. Isgur and M. B. Wise, Phys. Lett. B 232, 113 (1989); ibid. 237, 527 (1990). 



